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abstract 

We report on the discovery of Very High Energy (VHE) 7-ray emission from the BL Lacertae object 
lES 1011-1-496. The observation was triggered by an optical outburst in March 2007 and the source 
was observed with the MAGIC telescope from March to May 2007. Observing for 18.7 hr, we find 
an excess of 6.2<7 with an integrated flux above 200 GeV of (1.58±0.32) • 10~^^ photons cm~^ s~^. 
The VHE 7-ray flux is > 40% higher than in March- April 2006 (reported elsewhere) , indicating that 
the VHE emission state may be related to the optical emission state. 'We have also determined the 
redshift of lES lOlH-496 based on an optical spectrum that reveals the absorption lines of the host 
galaxy. The redshift of z = 0.212 makes lES 1011+496 the most distant source observed to emit VHE 
7-rays to date. 

Subject headings: gamma rays: observations, quasars: individual (lES 1011+496) 
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1. INTRODUCTION 

Known Very High Energy (VHE defined as > lOOGeV) 
7-ray emitting Active Galactic Nuclei (AGN) show vari- 
able flux in all wave bands. The relationship between the 
variability in different wave bands appears rather compli- 
cated. The MAGIC collaboration is performing Target of 
Opportunity (ToO) observations whenever alerted that 
sources are in a high flux state in the optical and/or X- 
ray bands. Previously, optically triggered observations 
resulted in the discovery of VHE 7-ray emission from 
Markarian 180 (| Albert et al.ll2006D . Here we report the 
discovery of VHE 7-rays from lES 1011+496 triggered 
by an optical outburst in March 2007. Previous obser- 
vations of the sourc e with the MAGIC t elescope did not 
show a clear signal (jAlbert et al.ll2007£l l. 

lES 1011+496 is a high frequency peaked BL Lac 
(HBL) object for which we now determined a redshift 
of z = 0.212 + 0.002 (Fig. 1). Previously, this was uncer- 
tain since it was base d on an assumed association with 
the cluster AbeU 950 (jWisniewski et al.|[l985,) . The red- 
shift determined here makes lES 1011+496 the most 
distant VHE 7-ray source yet detec ted with the possi- 
ble exception of PG 1553+113 (Ah aronian et al.ll2006al : 
[Albert et a l."2007b') (for which the redshift is 0.09 < z < 
0.42 (|Sbarufatti et al~l [20061 : iMazin fc Goebe]|[2?)0l L 

The spectral energy distribution (SED) of BL Lac ob- 
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Fig. 1. — Optical spectrum of lES 1011+496 obtained with 
the Multi Mirror Telescope, using the Blue Channel Spectrograph 
with the 300 line/mm grating, a 1.5" slit, and Loral 3kxlk CCD. 
Integration time was 30 minutes. Absorption lines of the host 
galaxy (Ca H&K, CH G, H^g and Mg b) are clearly visible and 
indicate a redshift of z = 0.212 ± 0.002. 

jects normally shows a two-bump structure. The lower 
frequency peak is due to synchrotron radiation. Various 
models have been proposed for the origin of the high- 
frequency peak; the most popular invoke inverse Comp- 
ton scattering of ambient soft photons. There have been 
several suggestions for the origin of the low-frequency 
seed photons that are up scattered to 7-ray energies: 
the soft photons may be produ ced within the j e t itse lf 
by synchrotron radiation (SSC, iMaraschi et al.l ()1992D ) 
or come from outs ide the je t, perhaps from the accre- 
tion disk (EC, DermeriLMilickeiscr (1993)). Thehigh- 
ener gy peak may, ins t ead, a lso have a hadronic origin 
(e.g. iMannheim et al.l (|199lD ). 

When the synchrotron emission peak is located in 
the low energy range from the sub-millimeter to opti- 
cal, the objects are called low-frequency-peaked BL Lac 
objects. HBLs, on the other hand, have the peak syn- 
chrotron emission in the UV to X-ray energy range. The 
peak of the second bump is often not observable be- 
cause of the low sensitivity above a few hundred MeV 
of satellite-borne detectors or a too high energy thresh- 
old of groun d-based 7-ray detectors. Wit h the excep- 
tion of M 87 (lAharonian et al.|[200l l2006b[ ) and BL Lac 
(| Albert et al.H2007d ). all known blazar sources detected 
at TeV energies with Cherenkov telescopes show a syn- 
chrotron peak in the UV to X-ray energy range, suggest- 
ing that the intensity of the TcV emission is related to a 
synchrotron component extending to high frequencies. 

2. OBSERVATIONS AND DATA ANALYSIS 

The MAGIC telescope is located on the Canary Island 
La Palma (2200 m above sea level, 28°45'N, 17°54'W). 
The accessible energy range spans from 50-60 GeV (trig- 
ger threshold at sma ll zenith angles) up to tens of TeV 
(| Albert et aLll2007dr i. 

The MAGIC observation was triggered by an observed 
high optical state of lES 1011+496 on 2007 March 12 
(see light curve Fig. 2). The source has been monitored 
for more than 4 years in the optical with the KVA^^ 
and Tuorla 1 m telescopes as a part of the Tuorla blazar 
monitoring program.^* In March 2007 the flux reached 
the highest level ever observed during the monitoring. 
The core flux, which is the host-g alaxy-subtract e d flux 
(the host galaxy flux is taken from lNilsson et all (j2007l ) 
and is 0.49 ± 0.02 mJy), increased more than 50% from 
the local minimum of the light curve. The high optical 
state with increasing flux was continuing throughout the 
MAGIC observations, despite an observation gap of 3 

http://tur3.tur.iac.es 
^® http://uscrs.utu.fi/~kani/lm/ 
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Fig. 2.— The optical R-band light curve of lES 1011+496 from 
Tuorla 1 meter and KVA 35cm telescopes. The vertical lines in- 
dicate the starting point of the MAGIC observations in 2006 and 
2007. 

weeks due to bad weather. 

lES 1011-^496 is monitored by RXTE ASM and Swift 
BAT, but the X-ray flux of the source is below the sen- 
sitivity of these instruments and the light curves show 
no indication of flaring. The source was also observed at 
Metsahovi Radio Observatory in May 2007. The source 
was not detected at 37 GHz, which indicates that it 
was not in a high state at millimeter wavelengths (A. 
Lahteenmaki 2007, private communication). 

After the alert, MAGIC observed lES 1011+496 in 
March-May 2007. The total observation time was 26.2 
hr, and the observation was performed at zenith angles 
ranging from 20° to 37°. T he observation was done in 
the so-called Wobble- mode ()Daum et al.l |1997[ ). After 
removing runs with unusual trigger rates, mostly caused 
by bad weather conditions, the effective observational 
time amounts to 18.7 hr. 

The data were analyzed using the standard anal- 
ysis a nd calibratio n programs for the MAGIC tele- 
scope (|Albert et al.i 2007d) . Th e analysis is based on 
image parameters (iHillad il985l) . the Random Forest 
(lBreiman|[200l iBock et all 120041) and the DISP meth- 
ods (jPomingo-Santamaria et al.l 120051 ) . After cuts for 
7/hadron separation, the distribution of the angle 9, 
which is the angular distance between the source posi- 
tion in the sky and the reconstructed shower origin, is 
used to determine the signal in the ON-source region. 
Three background (OFF) regions of the same size are 
chosen symmetrically to the ON-source region with re- 
spect to the camera center. The final cut 9^ < 0.02 deg^ 
to determine the significance (Fig. 3) was optimized on 
nearly contemporaneous Crab data to determine the sig- 
nificance of the signal and the number of excess events. 
The energy threshold was about 160 GeV for this analy- 
sis, which, given the soft spectrum of the source, allowed 
for signal extraction down to 100 GeV. The data were 
also analyzed with an independent analysis. Within the 
statistical errors the same significance, fiux, and differ- 
ential spectrum were obtained. 

3. RESULTS 

The distribution of the 0^-values, after all cuts, is 
shown in Fig. 3. The signal of 297 events over 1591 
normalized background events corresponds to an excess 
with s ignificance of 6.2cr using equation (17) inlLi fc Mai 
(fT98l . 

To search for time variability the sample was divided 
into 14 subsamples, one for each observing night. Fig. 4 
shows the integral flux for each night calculated for a 
photon flux above 200 GeV. The energy threshold has 
been chosen to reduce systematic effects arising from a 
rapidly decreasing effective area for 7-rays for lower ener- 
gies. The flux is statistically constant at an emission level 
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Fig. 3.— Distribution of 9^ for ON -region data and normalized 
OFF-rcgion data. The signal region is marked by the dashed line. 
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Fig. 4.— The night-by-night Ught curve of IBS 1011-1-496 from 
2007 March 17 (MJD 54176) to 2007 May 18 (MJD 54238). The 
year-by -year light curve is s hown in the inset, the 2006 data point 
is from [Albert et al.l ll2007al '). 



of F(>200 GeV) = (1.58±0.32) • 10"" photons cm-^g-i. 

The energy spectrum of lES lOlH-496 is shown in 
Fig. 5. It extends from ~ 120 GeV to ~ 750 GeV and 
can be well approximated by a power law: 



dA^ 
dE 



= (2.0 ±0.1) • 10" 
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The errors are statistical only. We estimate the sys- 
tematic uncertainity to be around 75% for the abso- 
lute flux level and 0.2 for the spectral index. The ob- 
served spectrum is affect ed by the e volvin g extragalactic 
background hght (EBL, iNikishovl (| 19621 ): IStecker et all 
([1992)) as the VHE 7-rays are partially altered by inter- 
actions with the low-energy photons of the EBL. There- 
fore, to obtain the intrinsic spectrum of the source, the 
observed spectrum must be corrected. The optical depth 
and the resulting attenuation of the VHE 7-rays from 
lES 1011-1-496 are calculated using the number density 
of the evolvin g EBL provided by the best-fit model of 
iKneiske et al] (|2002[ ). Within given model uncertainties, 
the model is in good agreement with alte rnative models 
(jPrimack etIaD|2005|; [Stecker et all 1 .20061) and EBL ud - 
per hmits (|Aharonian et al.i r2006c: M azin fc Rau3l2007f ). 
Even after the correction, the slope of the spectrum is 
Pint = 3.3 ±0.7 (dashed brown hue in Fig. 5, / NDF = 
2.55/2), softer than observed for other HBLs in TcV en- 
ergies and thus not providing new constraints on the EBL 
density. 



Fig. 5. — The measured spectrum (black filled circles), the power- 
law fit to the data (solid line), the deabsorbed spectrum (brown 
open circles), and the fit to the deabsorbed spectrum (dashed 
brown line). The last measured point is a 95% upper limit. In the 
deabsorbed spectrum, the last spectral point at rjOOOGcV is 1.6 cr 
above the fit and t hus not significant. T he Crab Nebula spectrum 
(green dashed line. [Albert et aTl Il2007dl) ) is shown for comparison. 

4. DISCUSSION 

We report the discovery of VHE 7-ray emission from 
BL Lac object lES 1011+496. With the redshift of z = 
0.212, it is the most distant source detected to emit VHE 
7-rays to date. The observed spectral properties (soft 
and no significant excess above ~ 800 Ge V) are consis- 
tent with the state-o f-the-art EBL models (Knciske et al.l 
120021 : iPrimack et"ai. 2005; Stc_ck cr.ct al, 2006) and con- 
firm recently derived EBL limits ( Aharonian et al.ll2006cl : 
iMazin fc Rau3 [20071) . 

In Fig. 6 we show the SEP of lES 1011±496 using 
historical data (open circles; ICostamante fc Ghisellinil 
(2002) and references therein) and our nearly simulta- 
neous optical i?-band data (triangle), together with the 
MAGIC spectrum corrected for attenuation (filled cir- 
cles). We also display (square) the EGRET flux of the 
source 3EG J 1009+4855, which has been suggested to be 
associated with lES 1011+496 (Har tman et al. (1999), 
but see also ISowards-Emmerd et alTI (|2003f ) whose anal- 
ysis disfavors the association). 

We model th e SEP by using a one -zone synchrotron- 
SSC model (see lTavecchio et al.l ()2001f ) for a description). 
In brief, a spherical emission region is assumed with ra- 
dius R, filled with a tangled magnetic field of mean inten- 
sity B. The relativistic electrons follow a broken-power- 
law energy distributions specified by the limits 7niin and 
7max and the break at 7b. Relativistic effects are taken 
into account by t he Poppler facto r (5. 

As discussed in iTavecchio et al.l (|1998l ) , if the position 
and the luminosity of the synchrotron and SSC peaks 
are known and an estimate of the minimum variability 
timescale is available, it is possible to uniquely constrain 
the model parameters. Unfortunately, we do not have 
all the required information to accomplish this. In par- 
ticular, we have fixed the synchrotron peak by requiring 
that it reproduces the optical flux and the historical X- 
ray spectrum and we assume the SSC peak to be close 
to the MAGIC threshold. These choices minimize the 
required emitted luminosity, since a lower SSC peak fre- 
quency would require a higher SSC luminosity. 

We present two models. The first (solid line), assum- 
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Fig. 6.— Spectral energy distribution of lES 1011+496. The 
two different fits are done by varying the minimum electron en- 
ergy 7niin (see text). The other fit parameters are: R(radius 
of sphere)=10^^ cm, (S(Doppler factor)=20, B(magnetic field) = 

0. 15 G, 7max (maximum electron Lorentz factor)= 2 ■ 10^ 71, (break 
electron Lorenz factor)= 5 • 10* the slopes of the electron distri- 
bution ni =2 and 712 = 5 before and after the break energy, 
respectively, as well as rio (normalization of the electron energy dis- 
tribution) = 2-10'* cm"''. The model is not intended for describing 
the radio data, which is assumed to origin from a larger emitting 
volume to avoid an intrinsic absorption. 

ing an electron distribution extending down to 7min = 

1, clearly overpredicts the MeV-GeV flux measured by 
EGRET. In the second case (dashed line), which fixes 
the low energy limit at 7min = 3 • 10^ (leading to a 
"narrowing" of both the synchrotron and SSC bump, see 

iKatarzvnsk i et aP 1)20061 )) the model is compatible with 
the reported EGRET flux. It is evident that simulta- 
neous GLAST -yiKGlC observations of this source will 
provide important constraints on the model parameters. 

In both cases, the energy output of the SSC compo- 
nent (reaching observed values of L ~ 10'^^ erg/s) dom- 
inates over the synchrotron luminosity, implying a rela- 
tively low magnetic field, B = 0.15 G. In that case the 
source would be strongly electron dominated, since the 
magnetic energy density would be several orders of mag- 
nitude below that of the relativistic electrons. A larger 



synchrotron flux (limited by the non-detection by BAT 
and ASM) could alleviate the problem. Simultaneous 
X-ray and VHE observations are mandatory to further 
investigate this issue. We also note the fit Doppler factor 
{5 — 20) is rather high and should be verified by Very 
Long Baseline Interferometric observations. The fitted 
parameters are similar to those derived for other TeV- 
emitting BL Lacs. We note however that adopting uioA- 
els where the jet has a velocity struc t ure (e.g. models 
bvlGeorganoDOulos k. Kazanas I ()2003D :lGhiseUini et al. I 
(|2005D ) would considerably reduce the required Doppler 
factors. 

lES 1011-1-496 was previously observed with the 
HEGRA telescope array, resulting in an upper limit of 
F(E>1 TeV)< 1.8 • lO'^^ photons cm'^s"! (Aharonian 
et al. 2004), which is well above the detected flux 
we found. The source was also observed by MAGIC, 
as part of a systematic scan of X-ray-bright HBLs, in 
March-April 2006. Being in a lower optical state (the 
core flux was ~ 50% lower than that in March-May 
2007), the observations showed a marginal signal with 
3.5ct significance corresponding to an integral flux of 
F(>180 GeV) = (1.26±0.4) ■ lO"" photons cm'^s'^, i.e. 
- 40% ()Albert et al.ll2007af) lower than the detected flux 
in March-May 2007 (see also the inset i n Fig. 4). A 
simila r trend was also found for BL Lac ([Albert et al.l 
l2007c( ). where the observations during a lower optical 
state failed to detect VHE 7-rays. This seems to indi- 
cate that there is a connection between the optical high 
state and the higher flux of VHE 7-ray emission at least 
in some sources. To further investigate this possibility, 
follow-up observations of the detected objects as well as 
further observations of other AGNs during high optical 
states are required. 
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